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Abstract. High magnetic fields in neutron stars, B ~ 
10 11 — 10 13 G, substantially modify the properties of atoms 
and their interaction with radiation. In particular, the 
photoionization cross section becomes anisotropic and po- 
larization dependent, being strongly reduced when the ra- 
diation is polarized perpendicular to the field. In a num- 
ber of previous works based on the adiabatic approxima- 
tion the conclusion was drawn that this transverse cross 
section vanishes for frequencies w smaller than the elec- 
tron cyclotron frequency uj c = eB/(m e c). In other works 
(which employed a different form of the interaction op- 
erator) appreciable finite values were obtained, ~ cro7 1 
near the photoionization threshold, where <To is the cross 
section without magnetic field, and 7 = 5/(2.35 x 10 9 G). 
Since it is the transverse cross section which determines 
the properties of radiation emitted from neutron star at- 
mospheres, an adequate interpretation of the neutron star 
thermal-like radiation requires a resolution of this contro- 
versy. 

In the present work we calculate the atomic wave func- 
tions for both discrete and continuum states by solving 
the coupled channel equations allowing the admixture be- 
tween different Landau levels, which provides much higher 
accuracy than the adiabatic approximation. This enables 
us to resolve the above contradiction in favour of the fi- 
nite transverse cross sections at u> < uj c . Moreover, for any 
form of the interaction operator the non-adiabatic correc- 
tions appear to be substantial for frequencies u> <; 0.3w c . 
The non-adiabatic treatment of the continuum includes 
coupling between closed and open channels, which leads 
to the autoionization of quasi-bound energy levels associ- 
ated with the electron cyclotron (Landau) excitations and 
gives rise to Beutler-Fano resonances of the photoioniza- 
tion cross section. We calculate the autoionization widths 



of these quasi-bound levels and compare them with the 
radiative widths. The correlation of the open channels is 
responsible for the modification of the cross section above 
the Landau thresholds. The results are important for in- 
vestigations of the radiation emergent from the surface 
layers of neutron stars. 

Key words: atomic processes: photoionization - atomic 
processes: autoionization - magnetic fields - stars: neutron 
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1. Introduction 

The proper interpretation of the recently discovered sur- 
face radiation of isolated pulsars (e.g., Ogelman 1995) re- 
quires knowledge of the elementary processes in magnetic 
neutron star atmospheres (Pavlov et al. 1995). Among 
such processes, ionization of atomic hydrogen is impor- 
tant both conceptually, due to the simplicity of the hy- 
drogen atom, and practically, because of the presumably 
strong gravitational stratification of the neutron star at- 
mospheres. 

Although the photo- and autoionization of atoms 
have been thoroughly investigated at the magnetic field 
strengths B ~ 10 5 - 10 9 G (e.g., Merani et al. 1995, and 
references therein), none of these works can be directly 
extended to the typical pulsar field strengths B 10 11 G. 
First, the atomic wave functions at the pulsar fields are 
nearly cylindrical, while at B = 10 9 G the ground state 
atom is still nearly spherical. Second, the Landau levels 
at the pulsar fields are displaced much more distantly, 
which alters the spectrum qualitatively. These qualita- 
tive features arise when the parameter 7 = ?kj c /(2Ry) = 
5/(2.35 x 10 9 G) exceeds unity. Here oj c = eB/(m e c) is 
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the electron cyclotron frequency, and Ry = m e 4 /(2?i 2 ) = 
13.6 eV is the Rydberg energy. 

The photoionization processes at 7 > 1 have been 
considered in a number of papers: first analytically by 
Hasegawa & Howard (1961) and Gncdin et al. (1974), and 
then both analytically and numerically by Schmitt et al. 
(1981), Wunner et al. (1983b), Mega et al. (1984), Miller 
& Neuhauser (1991), Potekhin & Pavlov (1993) (here- 
after Paper I), Bezchastnov & Potekhin (1994) (B&P), 
and Kopidakis ct al. (1996) (KVH). All these papers, ex- 
cept for the two latter ones, employed the adiabatic ap- 
proximation for the atomic wave functions. B&P and KVH 
studied the effect of atomic motion across the field on the 
photoionization process. Here we consider the particular 
case of atoms at rest. In this case the approach of KVH, 
who did not include electron cyclotron excitations, reduces 
back to the adiabatic one. There was no non-adiabatic 
treatment of the final state by B&P or by KVH. 

The photoionization cross sections obtained in the 
cited papers appeared to be strikingly different. Some 
authors (e.g., KVH) concluded that, for the non-moving 
atom, the cross section for photons polarized perpendicu- 
lar to the field vanishes in the most important frequency 
range lo < ui c , while others (e.g., Paper I) presented fi- 
nite transverse cross sections which, although suppressed 
by the strong magnetic field, are still significantly larger 
than, e.g., the Thomson cross section. Since the proper- 
ties of radiation emitted from magnetic atmospheres are 
mainly determined by the transversely polarized photons 
(Pavlov et al. 1995), the different values of transverse cross 
sections should then result in quite different spectra and 
angular distributions of the radiation. 

The first goal of the present paper is to resolve this 
contradiction, which is of crucial importance for model- 
ing the neutron star atmospheres. We shall prove that 
the discrepancy is caused solely by using the adiabatic 
approximation and is eliminated as soon as this approx- 
imation is abandoned. The zero values of the transverse 
cross section arise from using the so-called velocity form 
of the interaction potential, which (contrary to the alter- 
native length form) leads to inadequate results when used 
in the adiabatic approximation. The two representations 
yield the same nonzero transverse cross sections when the 
more accurate non-adiabatic approach is applied. 

Our second principal goal is to extend the consider- 
ation of the photoionization to the domains of magnetic 
fields and photon energies where the adiabatic approxi- 
mation does not provide sufficient accuracy. For instance, 
more than a dozen radio pulsars have magnetic fields in 
the range ~ 10 10 — 10 11 G; the corresponding dimension- 
less fields, 7 ~ 4 — 40, although they exceed unity, are 
not sufficiently high to neglect the non-adiabatic effects 
on the atomic properties. The non-adiabatic effects are 
also expected to be considerable for photon energies not 
very small as compared to the electron cyclotron energy, 
Tilo c ~ 1.2(B/10 n G) keV (Paper I). The spectral flux is 



maximal at these X-ray energies for effective temperatures 
~4x 10 6 (i?/10 11 G) K characteristic of hot polar caps of 
radio pulsars. Such radiation has been observed from a 
number of pulsars by ROSAT, and new observations are 
under way with ASCA. For lower temperatures, more typ- 
ical at the surfaces of cooling neutron stars, the radiation 
of such energies is still observable from nearby objects. 

Among the non-adiabatic effects, particularly interest- 
ing are some new qualitative features of photoionization, 
which arise from including non-adiabatic terms in the fi- 
nal (continuum) state as well as in the initial state. In 
this approach, due to coupling to the continuum of lower 
Landau levels, the autoionization of quasi-bound atomic 
states, enters naturally into the consideration, giving rise 
to the so-called Beutler-Fano resonances below the thresh- 
olds of the electron cyclotron excitations. Comparison of 
the autoionization widths of the quasi-bound levels with 
the radiative widths shows that autoionization is impor- 
tant at relatively low fields, B ^ 3 x 10 11 G. Our results 
demonstrate that the corresponding spectral features are 
expected to be observable in the thermal-like spectra of 
neutron stars with the next generation of X-ray satellites 
(particularly, AST RO~E) . In addition, above the electron 
cyclotron thresholds, different ionization channels are no 
longer independent of each other (as they were in the adi- 
abatic approximation). The correlation of the channels is 
responsible for considerable modification of the photoion- 
ization cross sections, which should be observable in the 
spectra and light curves of the EUV/X-ray radiation of 
neutron stars. 



2. Atomic wave functions 

2.1. Coupled channel equations 

The early studies of the hydrogen atom in strong mag- 
netic fields have been based on the adiabatic approxima- 
tion (e.g., Canuto & Ventura 1977). In this approximation 
the wave function of the relative motion ip srj (r) is factor- 
ized into a transverse part & ns (r±) and a longitudinal part 
gfP(z), where z is the relative coordinate along the mag- 
netic field, and r± = (x, y). The transverse part is just the 
Landau function which describes the transverse motion of 
a free electron in a magnetic field, n being the Landau 
quantum number, and s the negative of the ^-projection 
of the electron orbital momentum [n > 0, s > —n). The 
quantum number r\ = ±1 refers to the z-parity of the 
wave function. If the atom does not move as a whole 
across the field, which is the case of interest here, both 
X] and s are exact quantum numbers, whereas n can be 
considered as a "good quantum number" at 7 > 1. The 
longitudinal wave function gfP(z) obeys a one-dimensional 
Schrodingcr equation with a potential obtained by averag- 
ing the Coulomb potential over the transverse coordinates. 
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The total energy of an atom is the sum of the longitudinal 2.2. Bound states 



and transverse energies, E = E" 



E ns = [n+ (n + s)m e /m p ] Tiui c , 



EL, with 



(1) 



where the term (n+s)m e /m p accounts for the finite proton 
mass (Herold et al. 1981). 

In more recent studies the atomic wave functions have 
been calculated using various numerical approaches (e.g., 
Xi et al. 1992, and references therein). At 7 <; 1, the most 
convenient approach has been proposed by Simola & Vir- 
tamo (1978) and developed by Rdsner et al. (1984) and 
Potekhin (1994) (Paper II). It is based on the expansion 
of the wave function over the complete orthogonal set of 
the Landau functions $ ras in the xy plane. Since, for non- 
moving atoms, s is an exact quantum number, the cor- 
responding sum drops out, and we arrive at the reduced 
expansion 



^(r)= £ * n >.(r ± )g%(z) 



(2) 



where n m i n = for s > , and n m ; n = —s for s < . The 
functions g n v (z) can be found from the following system of 
coupled differential equations obtained by substituting the 
expansion (|^) into the Schrodinger equation, multiplying 
it by <&* s (rjJ and integrating over r±, 



tf_ d^ 
2/i dz 2 



V° n ,(z)g s n nz), 



(3) 



where fi is the reduced mass. The set of the effective one- 
dimensional potentials is determined as 

V° n ,(z) = (ns\ - e 2 /r\n's) ± = ($„ s | - e 2 /r|^ s ) = 

W.(OW,n'(0(£ + z 2 /2a 2 m )- 1/2 ^ , (4) 



a m \/2 

where /„„/ are the Laguerre functions (Sokolov & Ter- 
nov 1968), a m = ae7~ 1 ^ 2 is the magnetic length, and 
ob = h 2 /(m c e 2 ) is the Bohr radius. Neglecting the non- 
diagonal (n' ^ n) effective potentials decouples the system 
and brings us back to the adiabatic approximation. The 
choice of the maximum Landau number n max for trun- 
cating expansion (||) is dictated by the desired accuracy 
— the more terms one retains the higher accuracy of the 
wave function is provided. 

For any given set of the conserved quantities s and n, 
each term in expansion (^) is referred to as a channel, and 
the system (3) is the system of coupled channel equations. 
The whole set of n max — n min + 1 channels is separated 
into two groups. The first one includes open channels n = 
n m i n , . . . ,no — 1, for which E > E^ s . The second group 
embraces closed channels n — uq, . . . , n maX ) for which the 
opposite inequality, E < E^ s , holds. 



In the adiabatic approximation, the right hand side in 
(3) vanishes and the channels are uncoupled. The bound 
states have negative longitudinal energies, i. e., correspond 
to the closed channels. The discrete energy levels are enu- 
merated at specified s, and n by a longitudinal quantum 
number v = 0, 1, 2, . . .; the states having v = are known 
as tightly bound, while those with v > 1 are hydrogen- 
like. The longitudinal quantum number fully determines 
the z-parity of the state, rj = (— l) u , so that the quantum 
number rj becomes redundant and is usually omitted. 

In the non-adiabatic approach, the RHS in (3) cou- 
ples the state \nsv) into other channels, n' 7^ n, com- 
prising both bound states (closed channels) and con- 
tinuum states (open). The wave functions ip n sv( r ) — 
J2n' ®n's{r)^g™? v {z) and the energy levels E nsv can be 
designated by a leading term n' — n of the expansion (^]) , 
if 7 is large enough. Only the states with n = n m i n can 
be truly bound. In other words, the energies E nsv corre- 
spond to truly bound states only when all the channels are 
closed. Coupling of the closed channels shifts the level en- 
ergies from the adiabatic values and admixes higher Lan- 
dau orbitals to the bound state. 

A non-adiabatic computer code for calculating the 
bound state wave functions of the hydrogen atom mov- 
ing in a strong magnetic field has been described in Pa- 
per II. Here we apply this code (for the particular case of 
the non- moving atom) to obtain the initial state \i) of an 
atom subject to photoionization. 

2.3. Continuum 

The final atomic state |/) in the photoionization process 
lies in the continuum, Ef > E^ B) i. e., at least one asso- 
ciated channel is open. All previous publications on pho- 
toionization in very strong magnetic fields treat the final 
state adiabatically, i.e. (r\f) = $ Jls (rj_) g^\z). The sim- 
plest approach for the adiabatic continuum wave function 
is the Born approximation, used by Gnedin et al. (1974) 
for the non-moving atom and generalized by KVH to the 
case of motion across the field. It assumes the longitudinal 
wave function of the final state to be of the form 



gv)( z ) = exp(±ifc n z), 



(5) 



where k n — y/2\i{Ef — E^ s )/Ti is the electron wave num- 
ber, ± denotes the direction of motion of the outgoing elec- 
tron. A more accurate approach (e. g., Paper I) replaces 
exp(±ifc„z) by a function g^'(z) numerically determined 
from the corresponding uncoupled equation; this function 
turns into (^) asymptotically, at \z\ — > 00 (up to a loga- 
rithmic phase factor). 

In the non-adiabatic approach, the open (continuum) 
channel is coupled to closed channels and to other open 
channels (if the number of open channels hq — rt m i n > 1). 
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This coupling of open and closed channels causes, in par- 
ticular, autoionization of the quasi-bound states (states of 
positive energy which would be bound in the absence of 
coupling — see Section 2.4). The coupling of different open 
channels can be conveniently treated in terms of the reac- 
tance matrix R nn > (Seaton 1983). In the strong magnetic 
field case, the i?-matrix can be introduced as follows. 

System (3) of n max — n min + 1 coupled equations has 
(for a given parity 77) n max — n m j n + 1 linearly indepen- 
dent solutions (basis vectors g^ S7J = [gi^V 7 , ■ ■ ■ ,gn2^]) 
satisfying physically meaningful boundary conditions. 
It is convenient to enumerate the solutions similarly 
to the channels: n — n min , n min + l,...,n max (then 
gW s v g^ sv 5 nn i when the coupling switches off). 
Each of these solutions forms a wave function tpi^j' (r) = 
J2n> ®n>s(r±_) g^ sv (z) which comprises n - n min open 
channels and n max — uq + 1 closed channels. For unbound 
solutions (n = n m ; n , . . . ,n — 1) of a given z-parity, it is 
sufficient to construct real longitudinal wave functions for 
positive z, which satisfy the following asymptotic condi- 
tions at z — > +00: 

^(n,rcal) ^ ^ ^ ^ ^ ^ + R ^ ^ ^ ^ ^ 

for n' = n m j n , . . . , no — 1, and <7„™' r (z) — > for n' = 
ng, ... , n max . (Hereafter, in this and the following subsec- 
tions, we omit the indices s and n for brevity.) In Eq. (^|), 

4> n '{z) = k n >z + (k n >a B m c / ^y 1 \n{k n ,z) (7) 

is the z-dependent part of the phase. With the continuum 
wave functions normalized to unity in an interval of length 
L z , the i?-matrix satisfies the following symmetry relation, 

^n'^nn' — knRn'n- (8) 

The real basis of the wave functions g^' rca ^ (z) is con- 
venient for calculations, but it still is to be transformed 
into the basis of outgoing waves, appropriate to photoion- 
ization. This is done by analogy with the usual theory 
(Seaton 1983). For an electron outgoing in the positive z 
direction, the following asymptotic condition at z — > +00 
holds for the n-th solution: 

#K° ut) (z) _ exp[i0 n /(z)] + S\ n , exp[-i</>„,(z)], (9) 

where is the Hermitean conjugate scattering matrix. 
Now, if we compose a (n - n min ) x (n max - n min + 1) 
matrix function G(z) of the elements g^j\z) (with n the 
first and n' the second subscript), then we can obtain a set 
of the outgoing wave functions from the matrix equation 

G (out) = + ii? )-l G (rcal)_ ( 1Q ) 

The matrix is expressed in terms of the R matrix as 

= -{i + m)- 1 {i-\R). 



Note that the longitudinal wave functions satisfying 
the asymptotic condition (^|) should be multiplied by a 
common factor in order to ensure the correct normaliza- 
tion of the wave function t/>w '' out )(r) . It follows from the 
unitarity of the ^-matrix that this factor equals (2L 2 ) -1 / 2 , 
where L z is the z-extension of the periodicity volume of 
the final state. 



2.4- Autoionizing states 

The adiabatic and exact approaches treat in a fundamen- 
tally different way quasi-bound states associated with ex- 
cited Landau levels. Since the adiabatic approximation 
allows no coupling between the Landau orbitals, a sepa- 
rate set of bound states appears below each Landau level 
n. These states, however, lie in the n' < n continuum, 
and may therefore decay into the continuum via two pro- 
cesses. The first one is spontaneous emission of photons, 
which broadens these levels significantly in strong mag- 
netic fields. This process has been thoroughly studied by 
Wunner et al. (1983a). The second one is autoionization, 
which could not be accounted for in the earlier work based 
on the adiabatic approach. In the non-adiabatic treat- 
ment, due to the coupling, the electron in a quasi-bound 
state can escape to infinity via any of the open channels 
n' < n. 

The quasi-bound autoionizing states manifest them- 
selves as resonances of the continuum wave function. For 
weaker magnetic fields (7 < 1) such states have been stud- 
ied, e.g., by Friedrich & Chu (1983). Near a quasi-discrete 
level a resonance condition is satisfied, leading to a great 
amplification of the longitudinal coefficient g n (z). Thus, 
the shape of the wave function near the origin becomes 
similar to that calculated in the adiabatic approximation 
for the quasi-bound state. At large distances, where g n (z) 
decreases exponentially, the contribution of the orbitals 
n' < n dominates, which can be interpreted as electron 
leakage from the quasi-bound state. 

The general theory of autoionizing states has been de- 
scribed, e. g., by Friedrich (1991). Here we briefly discuss 
it for the case of the strongly magnetized hydrogen atom. 
Let the electron energy E be close to an energy E^ at 
which there would be a bound state in the closed channel 
n in the absence of channel coupling. To consider cou- 
pling of the channel n to an open channel n! < n, one 
can retain the two corresponding terms in expansion (|^) 
and two equations in the system (3), assuming there are 
no other quasi-discrete levels close to E%~, Let g^(z) and 
§n^{z) be the two linearly independent solutions for the 
uncoupled open channel (e. g., with asymptotic behaviour 
g°? (z) ~ sin (2) + <5 ad ] andg^f(z) ~ cos[(j) n >(z)+5 ad \ at 
z — > 00, where <5 a( j is the adiabatic phase shift, and <p n i(z) 
is defined by Eq. (0)), and g%l(z) be the solution for the 
uncoupled closed channel. Then, by analogy with the non- 
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magnetic case, the solution of the two coupled channel 
equations can be presented in the form 



g n ,(z) = cos<5 c g#(z) + sin<5 c g$(z) 



g nv {z) = -sm8 c — ■ lc . g nv (z) 



r a ,/2- 

nu.n' I 



(ii) 



(12) 



where S c is an additional asymptotic phase shift due to 
channel coupling, 



tan S c 



r a ,/2 

E — E„,, 



In Eqs. (O) and ©, 



(13) 



(14) 



is the (partial) autoionization width of the quasi-bound 
state \nv), 



EW= / K,n'(«) dz 



(15) 



is the coupling matrix element, and E nv is the resonance 
energy (slightly shifted from E^, — see Friedrich 1991). 
Note that Eqs. ( |l4] ) and ( |l5| ) can be derived directly with 
the help of the usual perturbation theory. 

Thus, it follows from Eqs. ([ll]) and (|l^) that coupling 
of the closed and open channels distorts the open chan- 
nel wave function and admixes the bound state \nv) to 
the continuum. The strength of the admixture is given by 
the square of the amplitude in front of the bound wave 
function gf^(z) in Eq. (p"2|); its dependence on energy is 
determined by the resonance function (Breit-Wigner pro- 
file) which coincides with the derivative of the phase shift 
S c with respect to energy, 



dSc 
dE 



v/2 



{E-E nv f + {T* ,/2f 



(16) 



The closed channel n is most strongly coupled with the 
open channel n' = n — 1. If we want to include coupling 
with all the open channels, the total autoionization width 

can be evaluated as = J2n~=n ■ ^nvn' • ^ s we snan 
see, the interference of the coupled states leads to the 
Beutler-Fano resonances in the energy dependence of the 
radiative transitions. 



3. Interaction with radiation 

3.1. Matrix element of the interaction 

The cross section for ionizing an atomic state \i) into a 
continuum state |/) due to interaction with a photon of 



frequency to, wave vector q, and polarization vector e, can 
be written as (e.g., KVH) 



J_ Ry 

8a 3 Tilu 




|(/|M|z)| 2 CTTh , 



(17) 



where a = e 2 /(he), and M is the dimensionless interaction 
operator. Its "velocity" representation involves the kinetic 
momentum operator tv = p + (e/2c)Bxr: 



where 



M, 



(t) 



a B exp(iqrj 



■ e ■ 7r - i(q x e) 2 



(18) 



(19) 



and 5M^> denotes corrections of order (m e /m p ). It will 
be shown that these corrections can be appreciable, when 
the velocity representation is used. Approximately one can 
write 

2 



m e i ( e 
m~ aB ~h~ e ' V ~ 2c X r 



(20) 



The first term in the square brackets in Eq. (|19|) represents 
the electron current or "velocity" term in the interaction 
potential, while the second term corresponds to the inter- 
action of the radiation magnetic field with the magnetic 
spin moment of the electron. We have neglected spin flip 
transitions, which are unimportant at B <, 10 13 G, ac- 
cording to Wunner et al. (1983b) and KVH, and fixed the 
electron spin antiparallel to the magnetic field. Besides, we 
have omitted the term eq (B&P) from the square brack- 
ets, assuming the transverse polarization. 

Using Eq. (||) for both atomic states, we obtain 



(f\M\i) = 53<n,/Kn,s / |M|n / 1 a < >x|n / 1 t)|| 1 



(21) 



where the subscripts _L and || denote the transverse (cf. 
Eq. (4)) and longitudinal matrix elements, respectively. 
The inner (transverse) matrix element in Eq. (|2l]) can be 
calculated analytically, using the well known properties 
of the electron Landau quantum states (e.g., Canuto & 
Ventura 1977), so that (f\M\i) reduces to the sum of one- 
dimensional quadratures. 

3.2. Transverse polarization: analytical consideration 

Using the commutation relations for the Hamiltonian, the 
matrix element (f\M\i) can be rewritten in the "length 
form" (Paper I): 



(f\M^ 



MM = M (r) +8M {r \ 



M, 



(r) _ 



iae exp(i</r) 



2m r uj 



e ■ r \ 1 — 



huj 

2m c c 2 



q ■ 7r 

m e uj 



(q x e) 



(22) 
(23) 



(24) 



(i 
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Here, as well as in Eq. (18), 5M denotes corrections ~ 
m e /m p . 

For exact atomic states \i) and |/), the two representa- 
tions are equivalent. However, the equivalence (22) breaks 
down in the adiabatic approximation. The most striking 
difference occurs for transitions within the ground Landau 
state when the photons are polarized perpendicular to the 
magnetic field. It was recently confirmed by KVH that in 
this case (0, Sf\M^\0, Si)±_ — identically, if one neglects 
the small corrections due to 5M^> . One actually finds an 
exact cancellation of the contributions arising from the 
velocity and the spin interaction terms in Eq. (p9|). At 
the same time, the representation (24) leads to a non-zero 
cross section!!. Some analytical estimates help to resolve 
this apparent contradiction (see also Appendix A of Paper 

I)- _ 

Since series (|2l| ) converges rapidly at 7 ^ 1, we can 

expect that it is sufficient to keep only the leading terms 
in it. In the velocity representation, the zero-order term 
(0,s/|M (7r) |0,Si)_L vanishes (KVH). Retaining the first- 
order terms, we obtain 



(i>/l <M/|M w l°= s iUIMii 

(0,/|{0 )S/ |M o W |l, Si )x|l,i)|| 
(0,f\(0,s f \SM^\0,s i ) ± \0,i) ]] . 



(25) 



Which of the three terms in Eq. (25) dominates, depends 
on the magnetic field strength and transition considered. 

In the length representation, the zero-order term dom- 
inates: 



r) 



z)«(0,/|(0, S/ |Af( r) |0, Sj ) ± |0,i)|| 



(26) 



Let us consider, for simplicity, the dipole approxima- 
tion, q — ► 0. It has been shown in Paper I (see also 
Sect. 4) that this approximation can be safely used at 
hu < am c c 2 ln7 ~ 300 In 7 Ry. For 9 = 0, Eq. (25) then 
yields 



(/|M« 



2i- 



.«o,/|M>n + cV57)*., 



- e_«l,/|0,i>|| + CV*)*. /ia< - 



(27) 



where e± = (e x ± ie y )/y/2 are the cyclic components of 
the vector e, and 



C = ^(o,/|o,i) 



(28) 



The terms with £ represent the leading contribution due 
to SM^K Analogously, from Eq. @ we obtain 



(f\M (r 



2iv^- 



UJ, 



y/sJS SfySi+1 



1 The spin term (q x e) 2 , omitted in Paper I, cannot be re- 
sponsible for this result, as was assumed by KVH. Indeed, the 
transverse cross section in Paper I remains non-zero even for 
photons directed along the field, in which case the spin term 
turns out to vanish identically. 



Si S Sf , Si -i (0,/|0,«;|| 



(29) 



Conditions under which Eqs. (27) and (29) give similar 
cross sections are obtained by equating their right-hand 
sides: 



(0,/|l,i) || «VS7f^-^) (0,/|0,0y, 



(l,f\0,i) 



m, 



(0,/|0,i 



(30) 



(31) 



where Sf — Sj + 1 or = Sf + 1, respectively. The 
term m c /m p comes from SM^ through Eqs. (27), @. 
Note that the terms m e /m p can give a substantial cor- 
rection (especially near the photoionization thresholds, 
huj ~ In 2 7 Ry) in very strong fields, when 7 ln~ 2 7 is not 
negligible in comparison with m p /m e . 

Our analytical estimates (Appendix A) and numerical 
results (Sect. 4) confirm that the approximate relation- 
ships d30|), ( |3l| ) are indeed satisfied if u> ^ lo c . Thus, in the 
adiabatic approximation (for n = n' = 0), the length rep- 
resentation of the interaction matrix element, employed 
in Paper I, enables one to calculate the leading contribu- 
tion to the transverse cross section which is missed when 
the velocity representation is used. On the other hand, 
it makes no difference which representation is used if the 
post-adiabatic corrections are included. 

3.3. Interference of open and closed channels: Beutler- 
Fano resonances 

As we have discussed in Sect. 2.4, coupling of a closed 
channel n to open channels, particularly to the channel 
n — 1, results in a resonance admixture of the quasi-bound 
state to the continuum. Interference of the two final cou- 
pled states leads to Beutler-Fano resonances in the pho- 
toabsorption spectrum (e. g., Friedrich 1991). According 
to Eqs. ( pr| ) and ( |T2| ) (for n' = n— 1), the transition matrix 
element, for the final electron energy close to the energy 
of the quasi-bound state, can be presented as 



(f\M\i) = Mi cos<5 c (1 - q tan<5 c ) , 

where 

_ M 2 U f 
q ~ Mi T^/2 



(32) 



(33) 



is the shape parameter of the Beutler-Fano resonance 
(Uf = Un,, T) = ^ — see Eqs. @ and 

TBI) 



Mi - (gf-i^Kn-hsflMlmsi)^^ (34) 
is the matrix element in the absence of coupling, 
M 2 = -(rj/2I7 / )(^ w Kn-l > a / |M|n i ai)x| SBt . m )|| 
+ (9ns f u f \(ns f \M\n i s i ) ± \g niSiVi ) ll (35) 
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is the coupling correction. The cross section is propor- 
tional to 



6 — 



|Mi| 2 cos 2 5 c tan<S c | 2 
\M x \ 2 F(e,q) , 



where 



is the Beutler-Fano function, 



e = - cot5 c = 2{E - E nSfVj .)/T) 



B=10 13 G 



(36) 



(37) 



(38) 



is the reduced energy. The function F(e, q) tends to unity 
in the very far wings of the resonance, at |e| ^> max(l, |g|). 
If q is real (e. g., in the dipole approximation), it turns 
to zero at e = —q and is maximal (F max = 1 + q 2 ) at 
e = l/q. In the limit of very weak coupling, Uf — ► 0, we 
have q oc UJ 1 — + oo, q 2 T^ remains finite, and the Beutler- 
Fano resonance turns into the delta-function (superim- 
posed onto the bound-free continuum) which describes the 
bound-bound transition in the absence of coupling. 

4. Results 

We present here example photoionization cross sections, 
calculated in accordance with the approach described 
above. In all cases the hydrogen atom is assumed to be ini- 
tially in its ground state. Basically we have employed the 
velocity representation of the interaction matrix element, 
including into M' 1 ' all the non-dipole and finite-mass cor- 
rections (Paper I), as well as the spin term (KVH). The 
results are compared with those obtained with the length 
representation, the adiabatic, and dipole approximations. 
Some details of the numerical techniques used are given 
in Appendix B. 

Our first principal result is that at ui <C uj c Eqs. ( j30| ) 
and ([H]) do hold with a high accuracy (the difference does 
not exceed 1-2% at Tiw ^ 30 Ry), thus confirming the 
considerations of Sect. 3.2 and Paper I. In the calculations 
we have also included higher terms of the expansions (|J) 
and (21); however, the channels n, n' > n/ + l proved to 
be unimportant. We have also checked the importance of 
the spin term, and found it to be negligible, by switch- 
ing it on and off numerically. On the contrary, the cor- 
rections due to the finite proton mass are appreciable. If 
the terms m c /m p on the right-hand side of Eqs. (|30|), ( |3l| ) 
were omitted, then the inaccuracy would increase to 10% 
at 7 = 1000. Clearly, this would introduce an inaccuracy 
of about 20% in the cross section. 

Figure |l| shows the cross section as a function of the 
photon energy hu for photoionization from the ground 
level, at a magnetic field strength B = 10 13 G (hu> c = 
8509 Ry); the cross sections plotted are for photons prop- 
agating along the magnetic field (9 = 0°) with the right 
circular polarization (er+), and propagating perpendicular 



o 



o 




— 



-2 — 



-3 — 



100 300 1000 3000 
Photon energy (Ry) 

Fig. 1. Total cross sections of the photoionization of 
the ground state H atom at the magnetic field strength 
B = 10 13 G. The curves are labelled by the symbols +, _L, 
and ||, corresponding respectively to o + (right circular polar- 
ization at 9 — 0°), cr± (polarization vector transverse to B at 
9 = 90°), and cry (polarization along B at 6 — 90°). Numeri- 
cal results (solid lines) are compared with the dipole approx- 
imation (short-dashed lines) and with the adiabatic results of 
Paper I (long-dashed lines) 



to the field (9 — 90°) with the linear polarizations paral- 
lel (cry) and perpendicular (erjj to B. Short-dashed lines 
correspond to the dipole approximation. We see that at 
9 = 90° the non-dipole corrections are unimportant for 
(Tii and <r± at Tioj S5 10 4 Ry. This confirms the validity of 
the transverse dipole approximation in this energy range 
and justifies negleting the spin term, as it depends on the 
transverse wave vector only. On the other hand, the inac- 
curacy of the longitudinal dipole approximation becomes 
perceptible at hu <; (2 — 3) x 10 3 Ry, which is seen from 
the deviation of the short-dashed line from the solid one 
for er + . 

For comparison, the results of Paper I (length repre- 
sentation) which do not include the non-adiabatic and 
spin terms are shown by long-dashed lines. The agree- 
ment with the adiabatic approach of Paper I, which in- 
volved the length form of the interaction matrix element, 
is fairly good at Tiui ^5 10 3 Ry. Note that the velocity form 
would lead to vanishing ct+ and u±_ (KVH). 
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Fig. 2. Same as Fig. ^ for B — 10 G. The symbol ' — ' corresponds to the left circular polarization at 9 = 0°. Numerical results 
for cr_ and ctm are plotted with the dot-dashed line, numerical results for a + and a± with the solid line. The long-dashed and 
short-dashed lines correspond to the adiabatic (Paper I) and dipole approximations, respectively, a Circular polarization (<r±) 
at 9 — 0°. b Linear polarization (<j±, ctu) at 9 — 90° 



The agreement with the adiabatic results becomes 
worse as lu approaches lu c from below, which is caused 
by the growing role of the "side" terms (closed channels 
n > 0) in Eq. (21). Figure]^ demonstrates the cross sec- 



tions at B = lfF* G (Huj c = 850.87 Ry). We see that the 
adiabatic approximation (long dash) may serve only as 
an order-of-magnitude estimate at to ~ u) c , whereas at 
uj & (0.2 — 0.3)w c the agreement is good again. For the 
magnetic field chosen, the channel n = 1 opens at the 
(threshold) energies hu> = 862.79, 863.25 and 863.72 Ry 
for Sf = — 1, and 1, respectively (these would be the 
n = 1 thresholds for the left, right and longitudinal polar- 
izations in the dipole approximation). Immediately above 
the thresholds, the adiabatic approximation for er_ and 
a± becomes sufficiently accurate again, as has been pre- 
dicted in Paper I. Photoionization from the ground state is 
strictly forbidden for the left polarization below the n = 1 
threshold. Above the threshold, the corresponding cross 
section er_ appears to be strongly underestimated in the 
dipole approximation (see the left panel of Fig. ||). The 
importance of the factor cxp(iqr) in this case is explained 
by approximate coincidence of the photon and electron 
wave numbers, as discussed in Paper I. The dipole approx- 
imation is inadequate also for <7|| (right panel of Fig. ||) 
above the n — 1 threshold because it misses the channel 



n = 1, Sf = — 1, which gives the main contribution just 
at these energies. For this polarization, the adiabatic ap- 
proximation (long-dashed lines) proves to be inaccurate 
as well, mainly because of the admixture of the (open) 
channels n = and n = 1. 

The solid and short-dashed curves for the non- 
adiabatic cross sections presented in Fig. || do not include 
values within ~ 30 Ry just below the n = 1 threshold. 
Within this gap complex resonance structures arise in the 
photoionization cross sections, too narrow to be resolved 
in the scale of Fig. |[ An example of such structures is 
shown in an expanded scale in Fig. ^. The resonances 
are associated with the quasi-bound autoionizing states 
|1, st, Vf) admixed to the continuum |0, s/). Their shape in 
the dipole approximation (short-dashed lines) is typical for 
the Beutler-Fano resonances (see Sect. 3.3) with a nega- 
tive q parameter: the peaks (at E — Ei SfUJ + {2q)^ 1 T\ SfVf ) 
are followed by troughs (at E = E\ SjVj . — qT\ SfVf /2). If 
\q\ 1 for a given resonance (which is fulfilled for the reso- 
nances shown in Fig. |) , then the FWHM of the peak coin- 
cides with the autoionization width T* a (see Sect. 2.4), 
and the distance between the peak maximum and the 
trough minimum is ~ \q\Tf Sfl , f /2 3> rf s v . For instance, 
for the main peak of a + (E), which is associated with 
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Fig. 3. Same as Fig. |2|in the energy range near the Landau threshold, a Circular polarization (<r±) at 9 = 0°. b Linear transverse 
polarization (a±) at 6 — 90°. Superimposed on the same diagram, the dotted curves show the "bound-bound" absorption profiles 
calculated without allowance for coupling of the closed and open channels; the shapes of these profiles are determined by the 
radiative deexcitation of the upper levels (see text for details) 



the autoionizing state 1 1 , 1,0), we have r^ 10 ~ 0.01 Ry 
and q ~ —200. Resonances of the same nature have 
been obtained previously for the case of lower magnetic 
fields, 7 < 1 (cf. Delande et al. 1991; O'Mahoni & Mota- 
Furtado 1991). For strong magnetic fields (and, conse- 
quently, higher photon energies), non-dipole corrections 
change the shape of the resonances so that the shape be- 
comes different for different polarizations and angles 9. In 
particular, the parameter q is no longer real, which re- 
sults in non-zero values of the trough minima shifted from 
their positions obtained in the dipole approximation. An 
additional effect entering when we go beyond the dipole 
approximation is the appearance of additional dipole- 
forbidden peaks overlapping with the dipole- allowed reso- 
nances. For instance, for er^, the peak at Tilo ~ 856.4 Ry, 
associated with the quasi-bound state |1, 1, 0), is preceded 
by another one, at Tllu ~ 854.5 Ry, attributable to the 
state 1 1 , 0, 0) , transitions to which are dipole-forbidden 
(due to a "transverse dipole" selection rule, see also Ven- 
tura et al. 1992). Similarly, dipole-forbidden transitions 
(i.e. opening of the n = 1, s = 1 continuum) are responsi- 
ble for the jump of a+ at 863.7 Ry. 

The second, much weaker and narrower resonance (its 
autoionization width is only ~ 1.5 x 10~ 4 Ry) in Fig. 
3a is associated with the hydrogen-like quasi-bound state 



1 1,1,2). It is not observed in Fig. 3b because of the or- 
ders of magnitude stronger background absorption due to 
the transition to the continuum state nf — 1, Sf = — 1, 
which is allowed at 9 ^ and whose threshold lies below 
this hydrogen- like quasi-bound level. In fact, there exist 
other Beutler-Fano resonances, related to the transitions 
to more excited quasi-bound states, but they are much 
too weak and narrow to play any role in the computed 
spectrum, so that we do not display them here. For the 
same reason, we also do not show the autoionizing reso- 
nances of (Tii . For the longitudinal polarization, only tran- 
sitions to the odd states are allowed at 9 = 90°. Therefore 
only hydrogen-like autoionizing states can contribute to 
(Tii . The corresponding resonances, however, are extremely 
weak and narrow (r a < 10~ 4 Ry). 

Figure 4 shows the cross sections at a weaker field, B = 
10 11 G. This field strength was not considered in Paper I 
because the adiabatic approximation may become too 
crude in this case. The present non-adiabatic treatment 
allows us to include this (and lower) field strength(s) into 
the consideration. Several Landau thresholds appear in 
the observationally relevant energy range in this case. Fig- 
ure 5 demonstrates the resonances associated with these 
thresholds. 
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Fig. 4. Total cross sections of the photoionization of the ground state H atom at the magnetic field strength B — 10 11 G. Solid 
and dash-dot lines correspond to the right and left circular polarizations, respectively, at the incidence angle 8 — 0°; short-dash 
and long-dash lines correspond to the transverse and longitudinal linear polarizations at 6 = 90° 



The cross section for the right circular polarization 
at 9 — 0° (solid lines in Figs. 4, 5) reveals a relatively 
broad resonance below each n-th Landau threshold asso- 
ciated with the tightly bound autoionizing states |n, 1,0) 
(their peaks lie at hui — 87.6, 173.0 and 258.4 Ry, and 
the autoionizing widths are T a = 0.014, 0.016, and 0.014 
Ry, for n = 1,2,3, respectively). A sequence of weaker 
and narrower resonances is further seen to converge to 
a corresponding threshold. These are related to the even 
hydrogen-like autoionizing states |n, 1,2), |n, 1, 4), . . . 

Analogous features for the left polarization (dot- 
dashed curves), associated with the states |n, — 1,0), 
|n, — 1,2), . . ., are seen for n > 2. The states |1,— l,v) 
are not coupled to the continuum (if the motion across 
the magnetic field is neglected) and do not contribute to 
the photoionization cross section. On the contrary, cou- 
pling of the states |n, — 1,0) for n > 2 is relatively strong 
(e. g., r a _ 10 = 0.014 and 0.016 Ry for n = 2 and 3, re- 
spectively), and the corresponding resonances dominate 
in Fig. 5b, c. The cross section a± for the transverse po- 
larization at 6 = 90° (short dashes) shows autoionization 
resonances associated with both \n,l,v) (for n > 1) and 
\n, —1, v) (for n > 2) because e± = (e+ + e_)/\/2 is com- 
posed of both circular polarizations. 

For the longitudinal polarization (long dashes), tran- 
sitions to odd states are only allowed at 6 = 90°. The 



strongest resonance (at Tilo — 90.0 Ry) in Fig. 5a corre- 
sponds to the state 1 1,0,1) {T\ m ~ 0.0008 Ry). Other 
odd-state resonances (barely seen in Fig. 5a) are consider- 
ably weaker. At higher Landau thresholds, all resonances 
of this type are weak and narrow. 

The Beutler-Fano resonances in Figs. 3 and 5 were cal- 
culated assuming that autoionization is the only channel 
for decay of the quasi-bound states, so that other mech- 
anisms which could lead to additional broadening of the 
resonances can be neglected. However, an excited state can 
always decay via spontaneous emission of a photon, and 
it is known (e. g., Wunner et al. 1983a) that the rate of 
the radiative decay can be very high in strong magnetic 
fields. The relative importance of the radiative and au- 
toionization decays is determined by the relation between 
the autoionization width T a and radiative width T r . If 
r a ^> r r , then most electrons excited to the quasi-bound 
state rapidly escape into the continuum, and the shape of 
the photoabsorption resonances is determined by the au- 
toionization. In the opposite case, spontaneous emission 
occurs faster than autoionization, so that absorption of 
radiation at resonance energies is not accompanied by pho- 
toionization, and the shape of the resonance is described 
by the Lorentz profile of the width T r . Thus, it is impor- 
tant to compare T a and T r for a given level in order to 
evaluate which of the two processes is more important. 




The (total) autoionization widths for a few levels 
| v) and |2, s, v) are plotted as a function of B in Fig. 6. 
We see that they decrease with increasing B, being quite 
different for different autoionizing states. To understand 
qualitatively the behaviour of the widths, consider, for 
instance, the dependence of the coupling matrix element 
Uisv.o, Eq. ([l5|), which couples the n = 1 quasi-bound 
states to the n = continuum, on the dimensionless mag- 
netic field 7. A characteristic length, ~ clm = »B7~ 1 / 2 , of 
the effective potential Viq(z) in Eq. (|l5|), which determines 
the limits of integration over z, is much smaller than a typ- 
ical size, ~ an (~ 03/1117 for v = 0, and ~ clb{v + l)/2 
for v > 0), of the bound wave function <?f^(;z). Therefore, 
in the integrand of Eq. ([l5|), we have ~ oij 1 for 

the even quasi-bound states, and ffi^(^) ~ fltii (aM/<2||) 
for the odd states. The continuum wave function g^iz) 
depends on the product kz ~ (27) 1 / 2 (z/a B ) ^3 \/2. Tak- 
ing into account that a characteristic magnitude of Vio(z) 
is ~ e 2 /aM ~ (e 2 /ctB)7 1 ^ 2 , we arrive at the follow- 
ing estimates for the autoionization widths at 7 ^> 1: 
Ff si , oc 7 _1 / 2 ln7 for the tightly bound states, oc 7 -1 / 2 
for the hydrogen-like even states, and oc 7~ 3 / 2 for the 
hydrogen-like odd states. Similar behaviour is observed 
for the autoionization widths of the states with n > 1. 

The radiative widths are mainly determined by the 
transitions \nsi>) — > |n — 1, s + 1, u); at 7^ 1 the widths 



depend only on the Landau quantum number n and the 
magnetic field, 

r; = , (39) 

where a is the fine structure constant. In fact, (|3^) co- 
incides with the well-know cyclotron width for cyclotron 
transitions of free electrons (e. g., Daugherty & Ventura 
1978). We see from Fig. 6 that the radiative width exceeds 
all the autoionization widths at B > 6 x 10 11 G for n = 1, 
and B <; 4 x 10 11 G for n = 2. This means that at very 
high magnetic fields the quasi-bound states are destroyed 
via radiative decay rather than via autoionization, and the 
shape of the resonances of the photoabsorption cross sec- 
tion is determined predominantly by the radiative broad- 
ening. At lower magnetic field, the shape of the resonances 
may be determined by different mechanisms for different 
quasi-bound levels. For instance, T a 3> r r for the lead- 
ing resonances in Fig. 5 (e.g. associated with the |1, 1, 0), 
1 2,— 1,0), 1 2, 1,0) states), whereas weaker resonances are 
subject to stronger radiative broadening. 

To illustrate the expected effect of radiative de- 
excitation on the photoabsorption spectra, we added to 
Fig. 3 Lorentz profiles of the "bound-bound" transitions 
to the states 1 1,1,0) and 1 1,-1,0) (dotted lines). The 
height and the shape of the profiles are determined by 
the radiative width (T\ = 0.19 Ry) and the oscillator 
strengths (/ooo,no = 3.1 x 10~ 5 , / O oo,i-io = 1.987). Note 
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Fig. 6. Autoionization widths of a few states |n, s, v) (labels near the curves) vs. magnetic field. For n = 2, the total widths, 
rLi/ = rfs^o + r2s„,i are shown. The autoionization widths of the states |2, —1, ^) coincide with those of the states |1, 1, v). 
The long-dashed curves show the radiative widths F r n of the levels with n = 1 and 2. The dots in the left panel show the values 
of the autoionization widths of the states (from top to bottom) |1, 0, 0), |1, 1, 0), |1, 0, 1) and |1, 1, 1) calculated by Friedrich and 
Chu (1983) for B = 2.35 x 10 9 G (7 = 1) 



that 1 1 , — 1 , 0) is the truly bound state (for atoms at rest), 
and transitions to it are actually the strongest amongst 
all the transitions from the ground state, so that this 
"cyclotron absorption" dominates in the photoabsorption 
spectrum near the n = 1 Landau threshold. It should 
also be mentioned that the transition |0,0,0) — > |1,1,0) 
is forbidden in the adiabatic approximation; its oscillator 
strength is provided by the admixture of 1 1,0,0) to the 
ground state, and |0, 1,0) to the excited state. 

5. Conclusions 

We have studied photoionization cross sections of the hy- 
drogen atom in magnetic fields B <~ 10 11 — 10 13 G, typical 
for pulsars. We have used exact interaction matrix ele- 
ments, including effects of finite proton mass, non-dipolc 
and spin interaction terms. Unlike previous authors, we 
use non- adiabatic wave functions for the initial and final 
states of the atom. This accurate treatment yields the fol- 
lowing conclusions. 

First, it resolves the acute contradiction of previous 
works concerning the cross section for photons polarized 
transversely to the magnetic field at energies Tiu < hu> c . 
This cross section is finite, being orders of magnitude 
larger than crxh near the threshold, in agreement with Pa- 
per I. This is shown to have no connection with neglecting 
the spin interaction, as was assumed by KVH. Moreover, 
at w < w c the present results nearly coincide quantita- 



tively with those in Paper I. We have proven that the zero 
values of a + and a±_ , obtained by a number of authors fol- 
lowing Schmitt et al. (1981), do not represent the reality 
but arise from their using the velocity representation of 
the interaction matrix elements in combination with the 
adiabatic approximation for the wave functions. This com- 
bination led those authors to miss the main contribution 
in the cross section. 

Second, the non-adiabatic treatment includes autoion- 
ization of the quasi-discrete levels associated with the 
Landau excitations. These levels are considered as truly 
discrete in the adiabatic approximation. Autoionization 
manifests itself in the photoionization cross sections as 
Bcutlcr-Fano resonances near the quasi-discrete levels; the 
shape of the resonances is determined by correlation be- 
tween closed and open channels. We have calculated the 
autoionization widths for the most important resonances 
in a wide range of magnetic fields and have shown that 
they exceed the radiative widths if the magnetic field is 
not too strong, B <, 10 11 G. 

Third, we have shown that the adiabatic results sig- 
nificantly deviate from the exact ones not only near the 
resonances, but in a rather wide range of photon energies, 
unless the condition w < w c is satisfied. Above the Lan- 
dau threshold, this deviation is due to correlation between 
different open photoionization channels. 

In this paper we restrict ourselves to the particular 
case of an atom which does not move across the field. It 
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has been shown by Pavlov & Meszaros (1993), B&P and 
KVH that transverse atomic motion leads to the opening 
of additional channels of ionization. In the particular case 
of slow motion, however, our present consideration allows 
one to decide between the small transverse cross sections 
presented by KVH and orders of magnitude larger cross 
sections obtained by B&P — in favour of the larger ones. 
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Appendix A: derivation of Eqs. ( |30[ ) 

In this section, we derive an approximate relation between 
non-adiabatic corrections which, according to Eq. (27), 
determine the interaction matrix element in the velocity 
representation for the transverse polarization, and the adi- 
abatic overlap integral (0, /|0, z)y which enters Eq. (29) for 
the matrix element in the length representation. 

Let us consider the quantum states \a) and \b) (n a 
rib = 0) belonging to adjacent s-manifolds: s a = s, 
Sb = s + 1. The longitudinal coefficients in Eq. (||) are 
given by the set of coupled differential equations (3) with 
the effective potential (4). At 7 ^ 1, the non-diagonal ef- 
fective potentials are small compared to the diagonal ones, 
and the non-adiabatic admixtures g n ^o are small respec- 
tively. Retaining only the leading terms in the equation 
for go(z), we have 



h 2 d 2 m 

77- -5-0 + E a - s huj c 

l[i dz z m p 

ft 2 d 2 „ , m ^ 

2[i az z m-n 



&\z) = V&Hz)gl> a) (z), (Al) 



(A2) 



Keeping the first-order terms in the equation for g± (z), 
we obtain 



s — - ) htjj, 

m r 



( a ) t \ 
9i (z) 



V^(z)g[ a \z) + V^\z)gi a \z). 



(A3) 



Here the terms V n ig n with n > 2 are neglected. 

Multiplying Eq. (Al) by g^ ] ' (z) and Eq. (A2) by 
g^\z), subtracting them by terms and integrating over 
z, we arrive at the approximate relation 



E„ — Eh 



(0,a\V Q 
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(A4) 



Analogously, from Eqs. (A2) and (A3) we obtain 



E a - E b - 1 



m e 
rn, 



(l,o|0,6)| 



(0,a|y ( 1 s) |0,6)|| + (l,a|^ 



0) 



v^ +1) I0,6)|| 



(A5) 



The second term on the right-hand side should be omit- 
ted in the approximation considered. Indeed, according to 
Eq. (4), the difference of two diagonal effective potentials 
can be expressed in terms of non-diagonal ones, 

Vtf(z)-Vti +I \z) = (s + l)- 1/2 

x yt+2V$(z)-V&\zj\, (A6) 

and thus the last term in Eq. (A5) represents a second 
order correction. Analogously, the difference on the right- 
hand side of Eq. (A4) equals 



V$(z) - V ^ +1 \z) = (s + l)- l < 2 V<Z\z). 

Comparison of Eqs. (A4) and (A5) yields now: 

(1, a|0, b) || — E a - E b + (m c /m p )hu) c 
~ Vs + 1 



= <0, a|0,6>]| 



E„ 



E b - (1 - m /m p )hu; c 



(A7) 



(A8) 



Applying this result to radiative transitions between a 
lower state \i) and an upper state |/), neglecting m e /m p 
in the denominator, and assuming lo <C w c , we arrive at 
Eqs. © and @. 

Appendix B: details of computation 

B.l. Solving the Schrodinger equation 

Bound state wave functions are calculated using the mul- 
ticonfigurational Hartree-Fock technique described in Pa- 
per II. 

Continuum wave functions are sought in the form 
(||), separately for each z-parity. The set of Eqs. (3) for 
the longitudinal coefficients is rearranged in two cou- 
pled subsystems, which are solved by two-step iterations. 
In the first step, the equations for the open channels 
(n, n' = n m in, . . . , no — 1) are solved, the contribution to 
the right-hand side from the other group of orbitals (n 1 = 
no, ... , n max ) being fixed. In the second step, the longitu- 
dinal wave functions of the closed channels (n, n' > no) are 
adjusted to the open-channel functions {n' < no) found in 
the preceding step. The procedure is then repeated; typi- 
cally, a few such iterations are sufficient to reach conver- 
gence. The only exception occurs in the resonance energy 
region, where the number of iterations increases up to 20- 
30, and the method finally fails in narrow energy regions 
corresponding to the very top of the peaks. Nevertheless, 
as can be seen in our figures, we are still able to trace 
substantial parts of the resonance profiles. 

The first step of each iteration is performed by the 
outward integration, employing the explicit fourth-order 
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Runge-Kutta scheme (e.g., Fletcher 1988) for the vector 
function g = (gj/\ . ■ ■ ,9^-i)- The integration extends 
to the point zq ~ 10 2 a m , where the off-diagonal effective 
potentials become negligible. The second step, solving the 
longitudinal equations for g n > no provided that g n <n are 
fixed, does not differ from that described in Paper II. After 
the iterative process ends, one extra integration for each 
orbital is required to proceed beyond Zo, where the orbitals 
are already uncoupled. 



B.2. Rearrangement 

The outward integration for the open channels is per- 
formed with the initial conditions arbitrarily chosen as 
<?n(0) = <?'j(0) = for all n < ng except n = j, whereas 
we chose gj(0) = 1 or g'j(0) = I depending on parity. In 
this way we obtain an arbitrary set of n$ linearly indepen- 
dent wave functions ^;v< arb ) (r), j = n m i n , . . . , no — 1. This 
set has to be rearranged, in order to meet the asymptotic 
conditions (||): 



9 



(n,rcal) 



(*)= E 



(j,arb) 
l nj9 n ' 



(*), 



(Bl) 



j=n u 



where the coefficients a n j constitute the rearrangement 
matrix A. In order to obtain A, first, the integration is 
extended to a point z as 3> e 2 / min(£'/ — ■ ), where the 
asymptotic behaviour is reached. The minimum longitu- 
dinal energy here is to be taken over all open channels, 
i.e. for n = 0, . . . , no — 1. An outermost part of this inte- 
gration may be performed by a faster scheme, cf. Paper I. 
Then, using Eq. ([Bl]) together with Eq. (g), we arrive at 
the algebraic system: 

no— l 

E -9n'' alb) ( z ns) a nj ~ COS (j) n / (z a s)-Rrm' 

J=ri m i n 

= 5 nn > sm</> n (z as ), 

= ^m'^( z as)cOS(/)„(2: as ), n, n' = 0, . . . , n - 1 , (B2) 

the phase <j) n being defined in Eq. (R). For each given n, 
the 2(n — n nlin ) x 2(n — n min ) system (B2) is solved 
to get the nth row of the matrices A and R. Since the 
matrix of this algebraic system does not depend on n, the 
FACT/SOLVE code (Fletcher 1988) is most useful. 

Longitudinal matrix elements, which enter Eq. (^l|), 
are calculated along with the functions <7<J' arb ) at the last 
stage of the integration for each specific j. They still need 
to be transformed into the matrix elements for the outgo- 
ing states. According to Eqs. (Bl) and (|Io|), this transfor- 
mation of the array of matrix elements is performed by the 
matrix \y/2/L z {l + ii?) _1 ^4, which acts on the subscript 
j related to a channel. The square root multiplier ensures 
the necessary normalization (Sect. 2.2). 
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